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ABSTRACT 

Context. Gamma-ray binaries could be compact pulsar wind nebulae formed when a young pulsar orbits a massive star. The pulsar 
wind is contained by the stellar wind of the O or Be companion, creating a relativistic comet-like structure accompanying the pulsar 
along its orbit. 

Aims. The X-ray and the very high energy (>100 GeV, VHE) gamma-ray emission from the binary LS 5039 are modulated on the 
orbital period of the system. Maximum and minimum flux occur at the conjunctions of the orbit, suggesting that the explanation is 
linked to the orbital geometry. The VHE modulation has been proposed to be due to the combined effect of Compton scattering and 
pair production on stellar photons, both of which depend on orbital phase. The X-ray modulation could be due to relativistic Doppler 
boosting in the comet tail where both the X-ray and VHE photons would be emitted. 

Methods. Relativistic aberrations change the seed stellar photon flux in the comoving frame so Doppler boosting affects synchrotron 
and inverse Compton emission difl"erently. The dependence with orbital phase of relativistic Doppler-boosted (isotropic) synchrotron 
and (anisotropic) inverse Compton emission is calculated, assuming that the flow is oriented radially away from the star (LS 5039) or 
tangentially to the orbit (LS I +6r303, PSR B 1259-63). 

Results. Doppler boosting of the synchrotron emission in LS 5039 produces a lightcurve whose shape corresponds to the X-ray 
modulation. The observations imply an outflow velocity of 0.15-0.33c consistent with the expected flow speed at the pulsar wind 
termination shock. In LS I +61°303, the calculated Doppler boosted emission peaks in phase with the observed VHE and X-ray 
maximum. 

Conclusions. Doppler boosting is not negligible in gamma-ray binaries, even for mildly relativistic speeds. The boosted modulation 
reproduces the X-ray modulation in LS 5039 and could also provide an explanation for the puzzling phasing of the VHE peak in 
LS 1-^61 = 303. 
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> 1. Introduction 



Gamma-ray binaries display non-thermal emission from radio to 
very high energy gamma rays (VHE, >100 GeV). Their spectral 
luminosities peak at energies greater than a MeV. At present, 
three s uch systems are known: P SR B 1259-63 (Aha ronian et al.l 
l2005bl). LS 5039 (lAharonianet al. 2005a) and LS I h-6 1=303 
jAlbert et al.ll2006l) . A fourt h system, HESS J0 632-H057 may 
also be a gamma-ray binary ("H inton et al.l 12009 ). The systems 
are composed of a O or Be massive star and a compact object, 
identified as a young radio pulsar in PSR B 1259-63. All gamma- 
ray binaries could harbour young pulsars (lDubusll2006l) . 

Electrons accelerated in the binary system upscatter UV pho- 
tons from the companion to gamma-ray energies. The Compton 
scattered radiation received by the observer is anisotropic be- 
cause the source of seed photons is the companion star. VHE 
gamma-rays will also produce e^e" pairs as they propagate 
through the dense radiation field, absorbing part of the pri- 
mary emission. This is also anisotropic. Both effects combine 
to produce an orbital modulation of the gamma-ray flux if the 
electrons are in a compact enough region. This modulation de- 
pends only on the geometry. Orbital modulations of the high- 
energy (HE, >100 MeV) and VHE fluxes have indeed been ob- 
served. The modulatioi is unambiguously identify the gamina-ray 
source with the bin ary dAharonian et al l2006llAlbert et alJ2006t 
lAcciari et al ] l2008h . 



Synchrotron emission can dominate over inverse Compton 
scattering at X-ray energies, providing additional information to 
disentangle geometrical eff'ects from intrinsic variations of the 
source. Suzaku and INTEGRAL observations of LS 5039 have 
revea l ed a stable modulation of the X-ray flux (Ta kahashi et alj 
12009^ iHoff^mann et al.l I2OO9I) . Possible interpretations are dis- 
cussed in §2. None are satisfying. The key point is that the X- 
ray flux is maximum and minimum at conjunctions and that this 
excludes any explanation unrelated to the system's geometry as 
seen by the observer. 



In the pulsar wind scenario, the synchrotron emission is 
expected to arise in shocked pulsar wind material coUimated 
by the stellar wind. This creates a cometary tail with a mildly 
relativistic bulk motion (Fig. [T]i. Relativistic Doppler boosting 
of the emission due to this bulk motion is calculated in §3 
with details given in Appendix A. The orbital motion leads 
to a modulation of the Doppler boost, as previously proposed 
in the context of black widow pulsars (Arons & Tavani 19931 
IHuang & Becked 120071) . The calculated synchrotron modula- 
tion is similar to that seen in X-rays in LS 5039. Although 
this is not formally confirmed due to their long orbital periods, 
LS I +61°303 and PSR B 1259-63 also appear to have modulated 
X-ray emission (Chernyakova et al. 2006, 2009; Acciari et al] 
I2OO9I lAnderfiub et al.ll2009h . The appUcation to these gamma- 
ray binaries is discussed in §4. 
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Fig. 1. Geometry of Doppler boosted emission from a collimated 
shock pulsar wind nebula. The orbit is that of LS 5039 (to scale). 
The comet tail moves away from the pulsar with a speed /3 - v/c 
at an angle ^aqw. If Onow - then intrinsic emission in the co- 
moving frame is boosted in the observer frame at inferior con- 
junction and deboosted at superior conjunction. 



2. The X-ray modulation in LS 5039 

2.1. X-ray observations 

LS 503 9 has shown steady, har d X-ray emis sion since i ts dis- 
coverv (iMotch et al.l IT9 97: Rib o et al.l Il999t iReig et al.; l2003t 
iMartocchia et al.ll2005HBosch-Ramon et a l."2005. 2007). RXTE 
obser vations hinted at orbital variabiUty (Bosch-Ramo n et aTl 
l2005h but co nfirmation had to wait the Suzaku and INTEGRAL 
observations dTakahashi et"ani2009HHoff^mann et al.ii200 9). The 
average spectrum seen by Suzaku from 0.6 keV to 70 keV is an 
absorbed power-law with spectral index a - 0.51 + 0.02 {Fy ~ 
V-") and A^H = 7.7 + 0.2 x 10^' cm^^ and Fi_io keV = 8 x 10"^^ 
erg cm"^ s"', consistent with previous observations. There is no 
evidence for a cutoff up to 70 keV. 

Variability in Suzaku is dominated by a well-resolved mod- 
ulation followed over an orbit and a half. The X-ray flux varies 
by a factor 2 with a minimum at (p 0.1, s lightly after superior 
conjunction (^sup - 0.05 based on lArago na et al. 2009) and a 
maximum at inferior conjunction (0inf - 0.67). The 1-10 keV 
photon index is also modulated, varying between 1.61+0.04 at 
minimum flux and 1.46+0.03 at maximum flux. The compari- 
son with Chandra and XMM measurements suggests the mod- 
ulation is stable on timescales of years (Kishishita et al. 2009). 
The column density is constant with orbital phase, as if there 
were only absorption from the ISM. The lack of significant wind 
absorption suggests that the X-ray source is located far from the 
system or that the wind is highly ionised and/or has a mass- 
loss rate g lO^^Moyr"' (Bosch-Ramon et al. 2007). Here, we 
assume that the X-ray source is situated within the orbital sys- 
tem. 




u (Hz) 




0.2 0.4 0.6 0.1 

orbital phase (periaslron=0) 



Fig. 2. Comparison of the model for LS 5039 described in 
§2.3 with observations. Top panel: spectral energy distribution 
showing the Suzaku 1- 10 keV maximum and minimum spectra 
dTakahashi et alj|2009[) flie 1 00 MeV - 10 GeV average Fermi 
spectrum ( Abdo et al. l2009bl) and the VHE spectra averaged 
over phas es INFC (dark points) and SUPC (grey points) as de- 
fined in Aharonia n et alj ([2006). The model spectra averaged 



over INFC and SUPC are shown as dark and grey curves. Middle 
and bottom panels: expected VHE gamma-ray and X-ray orbital 
modulation compared to the HESS and Suzaku observations. 



2.2. Inverse Compton X-ray emission? 

The phases of X-ray and VHE maximum (minimum) are iden- 
tical. If both are due to inverse Compton scattering off stellar 
photons then maximum emissivity is at superior conjunction. 
Subsequent in-system absorption due to pair production moves 
the observed VHE maximum flux to the inferior conjunction. 
X-ray photons are too weak for pair production but could be ab- 
sorbed in the stellar wind with a similar result. This can be ruled 
out since the modulation is seen in hard X-rays above 10 keV 
and A^H is constant with orbital phase. Thomson scattering of the 
hard X-rays is unlikely as it would require a column density of 
scattering electrons x 10"^ cm"^ (e.g. a Wolf-Rayet wind as in 
Cyg X-3 rather than an O star wind), two orders-of-magnitude 
above the observed absorbing column density and plausible stel- 
lar wind column densities. 
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2.3. Synchrotron X-ray emission? 

Alternatively, the X-ray emission is synchrotron ra diation from 
the same electrons that emit HE and VHE y-rays. In'Dubus et al." 
(2008), we proposed that several features of the VHE observa- 
tions could be explained by assuming continuous injection of a 
E^^ power-law of electrons at the location of the compact object 
in a zone with a h omogeneous magnetic field B of order 1 G 
dPubus e t al. "2008). The synchrotron X-ray spectrum expected 
in this modeQ is shown in Fig. |2l It is hard with a x 0.5. The 
electrons producing this X-ray synchrotron emission have ener- 
gies between 10 GeV and 1 TeV, for which the dominant cooling 
mechanism is inverse Compton scattering in the Klein-Nishina 
regime. This keeps the steady-state distribution close to the E'^ 
power law (Fig. 3 in Dubus et al. 2008). Synchrotron cooling 
takes over at higher energies, causing a break to a x 1 . In fact, 
the spectral index seen by INTEGRAL up to 200 keV is softer 
(a = 1 + 0.2) than the average index measured by Suzaku up to 
70keV (a = 0.51 +0.02). 

Whereas it is promising to have the hard X-ray spectral shape 
correctly reproduced, the level of X-ray emission is too low and, 
more importantly, the orbital X-ray lightcurve from the model 
is inconsistent with the observed modulation. The expected 1- 
10 keV lightcurve shows only a very modest change with a peak 
at periastron (Fig.|2]i. The reason is that the variations in particle 
and magnetic energy densities (a factor 4) compensate to keep 
the synchrotron emission almost constant. 

2.4. Variations in parameters? 

A better fit is possible by treating B or particle injection as 
free functions of orbital phas e or by taking adiabatic losses 
into account. iTakahashi et alj ( 12009 ) argued that the X-ray 
spectrum necessarily implies dominant adiabatic cooling of an 
E^^ electron distribution (this is sufficient but not necessary: 
as discussed above, Klein-Nishina cooling also keeps the 
distribution hard). The X-ray and VHE observations were then 
be fitted by adjusting the adiabatic timescale t^^ with orbital 
phase. The derived variation in fad mirrors the X-ray lightcurve 
with fad reaching a maximum at (p^. Ther e is no o bvious reason 
why fad should peak at this phase. Taka hashi et al .l (l2b09) expect 
the variation in fad to reflect variations in the size of the emitting 
zone, itself modulated by the external pressure of the wind. 
The relevant phases are those of apastron (low pressure) and 
periastron (high pressure), but not inferior conjunction which is 
an observer-dependent phase unrelated to wind pressure. In LS 
5039, 0inf is significantly diff'erent from the phases of periastron 
and apastron passage. Hence, it would require a coincidence for 
any intrinsic change in the source (B, number of particles, fad, 
size, etc) to result in a peak at this conjunction. 

The link between the extrema of the X-ray lightcurve and 
conjunctions calls for a geometrical explanation related to how 
the observer views the X-ray source. Doppler boosting (see 
Fig ID is a possible solution to this puzzle. 

' Here, the injected number of fresh particles is kept constant along 
the orbit whereas the energy density of cooled particles had been kept 
constant in Dubus et al. (2008). With the energy density constant, the 
particle distribution varied very little with orbital phase, which high- 
lighted the impact of anisotropic scattering. However, a constant in- 
jection in number of particles is probably more realistic (at least for a 
pulsar wind). It has no noticeable influence on the spectra but it slightl y 
changes the VHE lightcurve from that shown in iDubus et alj i2008l) . 
The VHE lightcurve remains compatible with the HESS results. 



3. Relativistic Doppler boosting 

In the interacting winds scenario, the X-ray emission is expected 
to occur beyond the shock where the ram pressures balance 
(Bignami et al. 1977; Maraschi & Treves 1981; Tavani et aTj 
1994; Dubus 2006). Particles in the shocked pulsar wind are ran- 
domized and accelerated. MHD jump conditions for a perpen- 
dicular shock and a low magnetisation pulsar wind give a post- 
shock flow speed of c/3 (Kennel & Coronitilll984t) . If the ratio 
of wind momenta rj - (Zi'p/c)/(M*v*) is small then the shocked 
pulsar wind is confined by the stellar wind. The shocked wind 
flows away from the companion star forming a comet-like tail 
of emission. Relativistic hydrodynamical calculations show the 
flow is conical with an opening angl e set by ?? and can reach 
highly relativistic speeds (Bogovalo v et alJl2008l) . High energy 
electrons emit VHE gamma-rays and synchrotron X-rays close 
to the pulsar and lose energy as they flow out, emitting in the ra- 
dio band far from the system (Dubus 2006). Here, the relativistic 
electrons radiating X-rays (by synchrotron) and VHE y-rays (by 
inverse Compton) are assumed to be localized at the compact ob- 
ject location. The calculation of the relativistic Doppler boosting 
in the flow is general and can also be applied e.g. to the case of 
a relativistic jet in a binary (Dubus et al...2010.) . 

3.1. Synchrotron 

Even if the flow is only mildly relativistic, Doppler boosting can 
introduce a geometry-dependent modulation of emission that is 
isotropic in the comoving frame (Fig. [TJ. This will be the case 
for synchrotron emission. The relativistic boost is given by 



r(l -ySCobs-eflow) 

where Caow is the unit vector along the direction of the flow and 
Cobs is the unit vector from the emission site, assumed to be the 
compact object location, to the observer. The flow will be as- 
sumed to be in the orbital plane where it makes an angle flflow to 
the star - compact object direction. 

The outgoing energy will be modified by e = Dobsf' and the 
outgoing flux will be Fyie) - D^^^^F[,{e'), with primed quantities 
referring to the comoving frame. In the case of a constant syn- 
chrotron power-law spectrum in the comoving frame with index 
a then 

^-syn - ^)^bT (2) 

The r atio of maximum to minimum flux is (see also lPelling et al.l 
II987h 

E^J\l^t\, (3) 

for /3-1/3, i-6Q°, q'=0.5. Relativistic boosting can significantly 
change the theoretical X-ray lightcurve discussed in §2. In the 
case of a purely radial flow (flgow^O"), maximum (minimum) 
boost occurs at the inferior (superior) conjunction ((/'obs - Jr/l-i 
or 7r/2 + i) where the flow is directed towards (away from) the 
observer. 

3.2. inverse Compton 

Inverse Compton emission will also be modified by relativistic 
aberration. The spectrum of the target photons seen in a given 
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solid angle in the comoving flow frame will be changed accord- 
ing to a different relativistic transform. If the star is assumed to 
be point-like, the relativistic boost involved is 



1 



r(l -jSe^.eflow) 
The total energy density from the star in the flow frame is 

2 OsbT^ 



47r 



(4) 



(5) 



with the Stefan-Boltzmann constant asb = 7.56 10 '"^ erg cm""* 
K The angle i// under which scattering occurs will also be 
changed. This angle (cos^' - ej^.e^,^^) is given in Appendix A. 
The inverse Compton spectrum is then calculated in the comov- 
ing frame as in Dubus et al. (2008). The resulting spectrum is 
then transformed back to the observer frame as in 33.11 

Because of this double transform, and because of the intrin- 
sic orbital phase dependence of scattering on stellar photons, the 
Doppler-boosted inverse Compton flux variability can be quite 
different from the Doppler-boosted synchrotron variability. In 
the case of Thomson scattering off a power-law of electrons 
dN cc -y^i'dy (see Appendix A) 



<x0;;+f(l-e*.eobs)' 



(6) 



Note that F^^ takes into account the decrease in target photon 
density with distance d to the star since the orbits are not circu- 
lar. Test calculations show this approximation captures the main 
features of the full calculation at hig h energies, including in the 
Klein-Nishina regime (see also Geo rganopoulos et al.ll200lh . It 
will be used to discuss the behaviour of the inverse Compton 
emission. 



4. Discussion 

The Doppler-boosted synchrotron (Fsy„, Eq. |2]l and inverse 
Compton (Fic, Eq. ^ intensity variations were calculated for 
the three gamma-ray binaries and are discussed here. Full cal- 
culations were also carried out for LS 5039 and LS I -h61°303 . 
The orbital parameters are t aken from [Manchester et a (11991 
for PSR B 1259-63 and from lAragona et all (l2009h forLS 5039 
and LS I +61°303. The inclination / is assumed to be /=30° for 
PSR B1259-63, and 60° for both LS 5039 and LS I -i-6r303 
(lDubusll2006h . 



4.1. LS 5039 

LS 5039 has a stellar wind velocity (v,,, 2500 km s"') signif- 
icantly greater than the compact object orbital velocity (vorb ^ 
400 km s"') so that the cometary flow is assumed to be purely 
radial (^aow - 0). Doppler boosting leads to peaks and troughs 
for the synchrotron emission F^yn at conjunctions as outlined in 
§3 (Fig.O. The amplitude of the inverse Compton flux (Fic) is 
reduced as the increased scattering rate at superior conjunction 
is compensated by a deboost of Dobs (and vice-versa at infe- 
rior conjunction). The shape of the modulation does not change 
much. The bottom panel shows that F^yn follows well the Suzaku 
data when /? is adjusted to 0. 15 in order to match the X-ray mod- 
ulation amplitude. The spectral index is fixed to the value ob- 
served by Suzaku, a - 0.5 (equivalent to p-2 for the electron 
distribution). However, this assumes the intrinsic synchrotron 
emission is constant with orbital phase, unlike what happens in 
the model discussed in §2 and shown in Fig.|2] 
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Fig. 3. Doppler-boosted synchrotron (Fsyn, Eq. |2]i and inverse 
Compton (Fic, Eq|6]l intensity variations for LS 5039 assuming 
/3=0 (top), yS = 1/3 (middle), y6=0.15 (bottom). In all panels a 
is 0.5 (equivalent to p-2) as given by the X-ray spectrum. The 
flow direction is radial (flflow - 0°). Dashed lines show Fi^ after 
attenuation due to pair production at 1 TeV. The bottom panel 
shows a comparis on of F^yg wi th the Suzaku X-ray measure- 
ments of Takahas hi et al ] (I2009h . The X-ray data is multiplied 
by a constant renormalization factor and y6=0.15 to match the 
X-ray amplitude. 



The precise relativistic corrections were applied to the model 
discussed in §2 (Fig.|2]l, assuming /3 - 1/3. No other changes 
were made. The average level of X-ray emission is not changed 
much. However, the relativistic corrections move the peak X-ray 
flux to superior conjunction and increase the amplitude of the 
variations, bringing the model X-ray lightcurve very close to the 
observations (shown in the bottom panel of Fig. |2]i. The spec- 
tral shape is slightly harder than the observed one by about 0. 15 
in the index a. The orbital modulation of a follows the X-ray 
lightcurve with a hardening of a from 0.42 (superior conjunc- 
tion) to 0.30 (inferior conjunction), which is similar in amplitude 
to the hardening observed by Suzaku (§2.1). However, the aver- 
age level of X-ray emission is systematically too low compared 
to the observations. Increasing the magnetic field by a factor 3 
would be sufficient to raise the level of X-ray flux but this would 
also modify the VHE spectrum, bringing the break at a few TeV 
to energies that are too low. The model assumes all the emis- 
sion arises within a single zone and this could explain this short- 
coming. The X-ray (and GeV) emission come from electrons 
that have significantly cooled since their injection and, thus, this 
emission would be more likely to be affected by a more detailed 
m odel wh ere particle cooling i s followed along the flow, as done 
in lDubusI (12006) based on the lKennel & Coronitil (11984 model 
for pulsar wind nebula. Numerical simulations are needed to pro- 
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Fig. 4. Same as Fig. 1 but with the corrections due to relativistic 
motion taken into account. The flow is assumed to originate at 
the compact object location with /3 - 1/3 and to point radially 
outwards from the star. 



vide detailed constraints on the geometry and physical condi- 
tions in the post-shock flow. 

As expected, the VHE gamma-ray lightcurve is not aff'ected 
much by the corrections because most of the escaping VHE 
gamma-rays are emitted close to inferior conjunction (as a re- 
sult of the 77 opacity). The modified VHE spectrum for SUPC 
phases is actually better than the original model that overesti- 
mated the VHE flux at a few TeV Pair cascading can fill in the 
flux between 30 GeV and a few TeV at this phase (Cerutti et 
al., submitted). At HE gamma-ray energies, in the Fermi range, 
the average flux level is reduced significantly because most of 
the HE gamma rays arise at superior conjunction where the 
flow deboosts the emission. Fermi observations of LS 5039 and 
LS I -i-61°303 show that the HE gamma-ray emission cuts off ex- 
ponentially at a few GeV, suggesting the emission in the Fermi 
range (100 M eV - 10 GeV) is a distinct component from the 
shocked flow (lAbdo et al.ll2009al b). This could be due to pulsar 
magnetospheric emission, in which case the relativistic correc- 
tions and model discussed here will not apply to the GeV com- 
ponent. 
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Fig. 5. Doppler-boosted synchrotron (Fsyn, Eq. |2]i and inverse 
Compton {Fic, Eq. |6]l intensity variations for LS I H-61°303 as- 
suming P-Q (top) and 1/3 (middle). Here, the direction of the 
flow is assumed to be tangent to the orbit (flflow 5^ 0). The bot- 
tom panel shows the inverse Compton emission above 100 GeV 
using the full calculation instead of Eq|6] Dashed lines show the 
inverse Compton emission after absorption due to pair produc- 
tion. In all panels a p-2 power-law of electrons, corresponding 
to a-Q.5 for synchrotron, is assumed. An off'set of 0.275 should 
be added to the above phases to compare with the radio-based 
ephemeris of LS I -H6r303. 

4.2. LS I +61° 303 

The impact of the relativistic Doppler corrections in 
LS I +6r303 (and PSR B1259-63) is more diflicult to as- 
sess because the orientation of the cometary flow is uncertain. 
The wind of the Be stellar companion is thought to be composed 
of a fast, tenuous polar wind and, more prominently, a slow, 
dense equatorial wind. These equatorial winds are effectively 
Keplerian discs with a small outflow velocity (compared to their 
angular velocity). If the compact object moves through this disc, 
then (neglecting corrections due to the orbital eccentricity) it is 
essentially moving through a static medium in the corotating 
frame, suggesting the outcome is more likely to be cometary 
flow trailing the orbit rather than directed radially away from 
the companion star. This wiU have to be confirmed by numerical 
simulations of the interaction. 

VHE observations by the MAGIC and VERITAS collabo- 
rations consistently find that the peak VHE emissi on occurs at 
phases 0.6-0.7 usin g the h istorical radio ephemeris tAcciari et al.l 
20081: lAlbert et al.lu009l) . The best estima tion of the periastro n 
passage phase in this ephemeris is 0.275 (lAragona et al.ll2009h . 
hence there is an offset of 0.275 between the radio ephemeris 
used by observers and the one used here. As outlined in §2, the 
phases of periastron/apastron passage or the conjunctions are the 
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natural phases where the physical conditions or the configuration 
of the system would be expected to produce minima or max- 
ima in the lightcurves. The peak VHE flux occurs 2 to 5 days 
before apastron and is clearly not associated with any of those 
phases, making it difficult to interpret only with anisotropic in- 
verse Compton emission and pair production. 

Superior conjunction in LS I -1-61° 303 occurs slightly be- 
fore periastron passage, and inferior conjunction slightly af- 
ter. The inverse Compton peak and trough match exactly 
with the conjunctions when there is no correction (top panel. 
Fig. |5]). Doppler corrections have a strong impact on the in- 
verse Compton lightcurve. Figure |5] shows the correction fac- 
tors for LS I -1-61° 303 if the flow velocity vector is taken to 
be exactly tangent to the orbit. The maximum boost is around 
phases 0.3-0.4 and the emission is deboosted around periastron 
passage. The effect is strong enough to push the maximum of 

and F^y^ at phases 0.57-0.67, using the radio ephemeris, 
as observed. The correlated behaviour is also consistent with 
the X -ray and VHE observations reported in lAnderhub et alj 
(120091). These conclusions also hold when doing a full calcu- 
lation (bottom panel. Fig. |5]) to properly take into account the 
Klein-Nishina cross-section. The calculation assumed a constant 
power-law distribution of electrons with p=2. The VHE spec- 
trum is Fy ~ v'^ because of Klein-Nishina effects and the X-ray 
spectrum is Fy oc v"*^^, both of which agree with observations. 

4.3. PSR B1 259-63 

The case of PSR B 1259-63 was also explored under the same 
assumption as LS I H-61°303 (Fig.|6]l. The inclination is relatively 
low / = 30° so that Fi^ is almost symmetric without Doppler 
corrections (top panel). Looking at the top two panels, it can 
be seen that the Doppler corrections have little impact on the 
overall lightcurve because of the low inclination. The bottom 
panel shows that high Doppler factors can strongly deboost the 
overall lightcurve even though the morphology remains roughly 
the same. There is no obvious relationship between these curves 
and the (sparse) X-ray or VHE observations. Other variability 
fac tors probably dominate in this much wider binary system. 

iBogovalov et al.l (l2008h carried out relativistic hydrodynam- 
ical simulations of a pulsar wind interacting with a stellar wind 
with the specific case of PSR B 1259-63 in mind. They found that 
the shocked pulsar wind can accelerate from bulk Lorentz factors 
X 1 close to the termination shock up to 100 far away. Emission 
from such highly relativistic flows is not compatible with ob- 
servations: the emission would be strongly deboosted (bottom 
panel, Fig.|6l) except where (and if) the line-of-sight crosses the 
relativistic beaming angle where it would produce a flare. The 
observed X-ray and VHE modulations in gamma-ray binaries 
suggest modest boosting. The X-ray and VHE emission is more 
likely to originate close to the termination shock where the jump 
conditions for an unmagn etized relativistic flow give /3 - 1/3 
(iKennel & Coronitilll984j) . 



5. Conclusion 

The X-ray orbital modulation of LS 5039 peaks and falls at con- 
junctions, suggesting that the underlying mechanism is related 
to the geometry seen by the observer. Phase-dependent Doppler 
boosting of emission from a mildly relativistic flow provides a 
viable explanation. The underlying assumption is that the flow 
direction changes with orbital phase, so that even constant in- 
trinsic emission becomes variable as seen by the observer The 




0.001 r 



0.5 1.0 1.5 2.0 2.5 

phase (periastron = 0) 

Fig. 6. Doppler-boosted synchrotron (Fsyn, Eq. |2]i and inverse 
Compton (Fie, Eq.|6]l intensity variations for PSR B 1259-63 as- 
suming /3-Q (top), 1/3 (middle) and 0.9 (bottom). The direction 
of the flow is assumed to be tangent to the orbit (flflow ^ 0). Note 
the logarithmic y-axis scale. Gamma-ray absorption is negligible 
here. 



peaks and troughs are at conjunctions for a flow directed radi- 
ally away from the star, as expected if the emission arises from a 
shocke d pulsar wind confined by the fast stellar wind of its com- 
panion (lDubusll2006h . A moderate relativistic speed of /3 - 0.15 
or 1/3 is enough to reproduce the morphology of the observed X- 
ray lightcurve a ssuming (resp.) eith er constant intrinsic emission 
or the model of iDubus et all (1200 8'). Note that these values of 
allow for quite large values of the opening angles. More detailed 
calculations assuming a conical geometry for the flow confirmed 
that the results were unchanged as long as the angular size of the 
flow is smaller than l/T (if larger, the modulation is dampened). 
Reproducing the level of X-ray emission is difficult with a one- 
zone model as it requires values of the magnetic field that are a 
factor 3 above current values, leading to cutoff in the VHE specta 
at energies that are too low. A more complex multi-zone model 
of the post-shock flow might resolve this discrepancy. 

Inverse Compton scattering in the flow of external stellar 
photons will be modulated differently than intrinsic emission 
from the flow. In the case of a radial outflow, the external seed 
photon flux will be deboosted at all phases. However, a flow tan- 
gent to an eccentric orbit, as might arise in LS I -f61°303 and 
PSR B 1259-63, can lead both to boosts and deboosts in the co- 
moving frame depending on orbital phase and thus give rise to 
complex modulations. The calculated Doppler corrected emis- 
sion in LS I +61°303 peaks in phase with the observed VHE 
maximum. This is noteworthy since a simple explanation had not 
yet been proposed for the phase of VHE (and X-ray) maximum 
in LS I -i-61°303. This explanation requires that the shocked pul- 
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sar wind flows along the orbit, which appears co mpatible with 
the rad io VLBI images on larger scales shown in lDhawan et alJ 
120061) . 

The present work assumed a pulsar relativistic wind in the 
orbital plane but microquasar models have also been proposed 
for both LS 5039 and LS I +6r303. In this case, the emission 
arises from a relativistic jet. The jet angle to the observer re- 
mains constant along the orbit and so do Dohs and F^yn- Hence, 
no orbital modulation of intrinsic (synchrotron) X-ray emission 
due to Doppler boosting would be expected, apart from the pos- 
sible impact of jet precession on timescale s longer than the or- 
bital period (Kaufman Bern ado et alj2002h . Doppler boosting in 
a relativistic jet cannot explain the X-ray modulation in LS 5039 
or LS I -i-61°303. However, unless the electrons are far from 
the system or the system is seen pole-on, the angle of interac- 
tion between photons and electrons e*.eobs will change with or- 
bital phase. A modulation in Fi^ is unavoidable. This variation 
in inverse Compton emission can explain the orbital modulation 
seen in high-energy gamma-rays from the microquasar Cygnus 
X-3 by F ermi Gamma-ray Space Telescope dAbdo et al. 1 120093 
iDubuse t al. 2010). 
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Fig. A.l. Geometry of the binary -I- pulsar wind nebula flow sys- 
tem. The calculations assume that the massive star is point-like 
and that emission in the tail is limited to a small region at the 
pulsar location. 
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Appendix A: Doppler boosted inverse Compton 
emission on stellar photons 

The star is approximated as a point source of photons and the 
electrons are confined in a very small region. The overall ge- 
ometry and vectors are shown in Fig. lA.ll In the point-like and 
mono-energetic approximation, the stellar photon density in the 
observer frame is 



dn 
dedO. 



(A.l) 



where eo is the incoming photon energy and //o is the cosine 
of the angle between the incoming photon and the electron di- 
rection. Applying relativistic transforms to go to the comoving 
frame gives 



dn' 



J = r^(l -ySe^.Cflow)- 



dn 



dn 



de'dO.' " dedQ. " dedO. 

Developing the Dirac functions leads to 
dn' 



de'dil' 



: n'fy6(e' - e^)6(fi' - h'q) 



(A.2) 



(A.3) 



with Hq = SD^^no and - D^'eo- For inverse Compton scat- 
tering on an isotropic distribution of electrons in the comoving 
frame, fi'^ w ^★■^obs < Dubus e t al. 2008,) . The unit vector trans- 
forms in the comoving frame as 



+ [(F - IXc.Cflow) - Fyg] Cflo 
F(l -y6e*.eflow) 



(A.4) 



The transform giving e'^^^^ is simply given by replacing e* with 
Cobs above. The dot product of the two vectors in the comoving 
frame simplifies to 



1 - e;.e't„ = £)obs£>*(l - e^.Cobs) 



(A.5) 



The anisot r opic inverse Compton scattering kernel in 
iDubus et~an ([2008) can then be used, with the photon den- 
sity given in Eq. IA.3I and with the direction given by e^.e^j^j,. 
The resulting outgoing spectrum is then transformed back to the 
observer frame by using ei = DobsfJ and Fy{ei) - £)|^|^^F^(e[) as 
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discussed in §3.1. £)obs is defined in Eq.[T|and ei is the outgoing 
photon energy. 

For inverse Compton emission by a power-law distribution 
of electrons in the Thomson regime, the spectrum in the comov- 
ing frame is given by 



F'M)^Kn',{l-eM^\'-}\ (A.6) 



where p is the power-law index and K is a constant. In this case, 
the spectrum seen by the observer is 

FM) = KnoD-:-"Dl^:(l - e^.e^b,)""' ( (A.7) 



so that, using the dot product in Eq. |A 

FAei) = &o2)o;f(l - e^.Cobs)"^' ( (A.8) 

where a = (p - l)/2 . This is identical to the expression fo und by 
iDermer et al. and ' Dermer & Schlickeised (Il993h in the 

case of external scattering by a jet propagating away from the 
seed photon source (an accretion disc). The formula in Eq. lA.61 
lA. Bl are formally only valid for Thomson scattering on an infinite 
power-law of electrons. 

For completeness, the orbital separation d is given by 

^-T^^^^. (A.9) 
1 + e cos(0 - ll)) 

with the semi-major axis a - (GMP^^^-^l'XTip-yi^, M the total 
mass, e the eccentricity, 9 the true anomaly and w the periastron 
angle of the compact object. If the flow is in the orbital plane 
where it makes an angle flflow to the star - pulsar direction then 

Cobs-Cflow = - sin(6i + flflow) sin / (A. 10) 

e^.Cobs = -sinflsin; (A. 11) 

where / is the incHnation of the system. 



